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Quenched Staggered Spectrum at j3 =6.0, 6.2 and 6.4 * 

Gregory Kilcup 3 
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We give a preliminary report on the hadron spectrum on an ensemble of quenched lattices at f3 values of 6.0, 
6.2 and 6.4, using staggered fermions and collaborators Rajan Gupta and Steve Sharpe. Because of the relatively 
small number of configurations we pay marginally more attention to the question of statistics than in previous 
analyses. We are unable to discredit quenched QCD. 
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1. The Ensemble 

As part of our ongoing program in weak matrix 
elements, we have accumulated a set of quenched 
gauge configurations and staggered quark prop- 
agators at several different (3 values. As a first 
step in our analysis we have computed the hadron 
spectrum. For the results on the scaling of matrix 
elements, see elsewhere in these proceedings. 

The parameters of the simulations are listed 
in Table 1. The configurations were generated 
principally using the overrelaxed algorithm, with 
one 20-hit Metropolis step out of every 10 sweeps. 
Staggered quark propagators were computed us- 
ing periodic boundaries in space, and both pe- 
riodic and antiperiodic in time. Taking linear 
combinations, we formed forward- and backward- 
going propagators. These propagators are exactly 
the same as would result from doubling the lat- 
tices in the time dimension, but this method saves 
memory space. For each quark mass and bound- 
ary condition we computed four propagators, us- 
ing two types of wall sources (with phases 1 and 
(— ) n ) and two source timeslices (0 and N t /2). 
The wall source was defined by fixing to Landau 
gauge, as this was convenient for matrix element 
calculations. 

2. Curve Fitting 

Having measured the average value of some cor- 
relator G(t), and chosen some fit function / (a; t), 
the generic curve fitting problems are (1) to de- 
termine parameters a which yield the best fit, 
and (2) to decide if the fit is reasonable or not. 



Table 1 

The Statistical Ensemble 



Beta: 


6.0 


6.2 


6.4 


Volume: 


16 3 x 24 


32 3 x 48 


32 3 x 48 


^samp ■ 


41 


23 


24 


^quark ■ 


.03 


.025 


.015 




.02 


.015 


.010 




.01 


.010 


.005 






.005 





The standard method in use in the lattice com- 
munity[l] is to solve both problems at once by 
estimating the covariance matrix from the data 
itself, and then minimizing the usual goodness- 
of-fit statistic 



(1) 
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Here SG is the vector of deviations between the 
data and the model curve, while C'(t, t') is the es- 
timated covariance matrix. In the limit of large 
sample size, the fluctations in the average value 
of a correlator converge to a Gaussian distribu- 
tion described by the covariance matrix, and the 
statistic £ is distributed as a \ 2 - One can then 
test the hypothesis that f(a;t) is the true infi- 
nite sample mean of value of G(t) by asking how 
often a random Gaussian fluctuation away from 
this hypothetical true mean would lead to a value 
of £ as large or larger than the one observed. This 
fraction is the confidence level of the fit, which is 
readily computable with standard routines. Typ- 
ically one demands a confidence level of at least 
.05 to trust a fit at all. 

For finite sample size TV, one may proceed by 
assuming that the fluctations SG about the true 
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Figure 1. Modes of the covariance matrix for the 
7r channel at [3 = 6.4 and m q = .015. Plusses and 
crosses indicate positive and negative values. 



mean are at least approximately Gaussian, and 
that ^ is approximately a \ 2 statistic. More pre- 
cisely, if we fit to T timelices and have a sam- 
ple of TV gauge configurations with N > T, then 
is distributed as an _F-statistic with T and 
N—T degrees of freedom[2]. (If x\ and are dis- 
tributed as x 2 with ri\ and rii degrees of freedom 
respectively, then by definition y = ^- follows 
an F-distribution with rii and rii degrees of free- 
dom.) For large N this of course reduces to an 
ordinary \ 2 5 but for of order T (which it always 
is for us), the F-distribution has a long tail. For 
samples small enough that T > N (as is too often 
the case), the covariance matrix is singular, and 
one typically proceeds by defining £ in terms of 
the singular value decomposition (SVD) inverse. 
In either case, treating £ as an ordinary \ 2 leads 
to "peculiar" fits where the fit curve misses all 
the data points, and gives misleading confidence 
levels. 

The root of the problem is that with N sam- 
ples, one cannot really hope to reliably determine 



N eigenmodes of the covariance matrix. These 
eigenmodes are notoriously ill-determined. In fig- 
ure 1 we show the first 11 modes of the covariance 
matrix for one of our pion correlators, scaling each 
unit-normalized eigenvector by the corresponding 
eigenvalue. As expected, the dominant mode of 
fluctuation is the one where all timeslices oscillate 
up and down together. The error bars are deter- 
mined by jackknife, and are drawn on only the 
first two modes, since beyond that they become 
huge. The most ill-determined modes have eigen- 
values several orders of magnitude smaller than 
the first mode, and these then dominate the fit. 
In figure 2 we resolve a typical uncorrelated fit 
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Figure 2. Resolution of a two particle fit in the 
space of (unit-normalized) eigenmodes. Plusses 
and crosses indicate positive and negative values 
of coefficients for the data; squares and diamonds 
show the fit function. 



in the space of eigenmodes. We see that most all 
of the lower modes are well fit, while the higher 
modes contribute hundreds of units of "x 2 " to the 
£ statistic. In sum, 
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• Minimizing £ can lead to pathological fits 

• Even a "good" fit may have artificially large 

x 2 

There are several cures to the problem of an 
unstable covariance matrix in small samples. Two 
unsatisfactory ones are: 

• Increase the statistics several-fold 

• Fit to points in a smaller time range 
Two more promising choices are 

• Fit the covariance matrix to a parametric 
form, e.g. the free field form for the appro- 
priate four-point function. 

• Fit to a smaller number of eigenmodes, 
which are reasonably well determined 

We have opted for the latter method, truncating 
the observed covariance matrix to a small (typi- 
cally 6-10) number of modes before inverting and 
defining a statistic 

?(M) = ^f Z^tie.-SGf (2) 

where TV is the sample size and M is the num- 
ber of modes kept. For M = T this of course 
is the usual F-statistic. One can verify this does 
something reasonable in a simple Gaussian model 
of our data. We construct a multi-normal distri- 
bution with zero mean and with covariance taken 
from typical real data. Drawing samples of N vec- 
tors from this sample, we then ask for the £ of the 
point 0, i.e. the true mean. In figure 3, we plot 
the confidence levels for such fits, i.e. the frac- 
tion of samples which have worse £ as a function 
of the reduced £. For the parameters shown in 
the figure, the conventional SVD goodness-of-fit 
statistic is large, and if interpreted as a x 2 would 
lead to the overwhelming rejection of the hypoth- 
esis that is the true mean at most any confidence 
level. By contrast, the mode-truncated statistic 
appears to be close to a x 2 - In practice we take 
advantage of this and assume it is a \ 2 when as- 
signing confidence levels our fits. Of course it 
would be better to generate the correct distribu- 
tion, but the difference is not crucial. 
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Figure 3. Comparison of the conventional SVD 
X 2 statisitic and the reduced mode statistic, using 
a realistic spectrum of eigenvalues. 



As a practical matter, one must choose how 
many modes to keep by studying the effect on the 
bottom-line fit parameters. As shown in figure 4, 
the dependence is typically very weak until too 
many modes are added and the confidence level 
drops to near zero, a which point we reject the 
fit anyway. Based on such studies we typically 
choose 6-10 modes depending on the channel and 
fit range. 

Secondly one must choose a fit model. We 
would have preferred to fit to the entire time do- 
main, using several particles in each channel, or 
including some perturbative model of the short 
time behavior of our correlation functions. How- 
ever, because we use extended (wall) sources we 
cannot claim to know the short-time behavior. 
Further, since we use Landau gauge, our wall 
sources evidently create correlations which can in 
part be understood at the independent propaga- 
tion of quark states. To model this early time be- 
havior properly would require too many param- 
eters. By using simulated annealing and other 
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minimization techniques, our fitting routines can 
handle models with a few (e.g. three) poles, but 
even this is not enough. Accordingly we chose to 
stick to single particle fits, and used the following 
protocol: 

• For each time t m ; n perform fits the region 
^min-^iast varying A^ mo de between 6 and 12. 

• Find the smallest t m ; n for which a "reason- 
able" (CL > 0.2) fit can be found. 

• Throw away one timeslice and one mode for 
the definitive fit. 

This fitting was performed automatically, and 
statistical errors were assessed using 250 boot- 
strap samples. In every case quoted, the errors 
were essentially symmetric, so we give only a sin- 
gle number for the error. 
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Figure 4. Mass parameters in a two particle fit 
to the 7r correlator as a function of the number of 
modes kept. Fits to the right of the dashed line 
have unacceptable confidence levels. 



3. Results 

The main results are summarized in Table 2, 
where we show various parameters extrapolated 
to the chiral limit. The first two rows should be 
compared to the physical values of 1.22 and 10.1 
respectively. Evidently the bottom line does not 
look too bad at f3 = 6.4. The parameter A is the 
mass of the the non-Goldstone pions in the chiral 
limit, obtained by fitting m\ = A 2 + Am q . The 
ratio A/m 2 p remains fairly constant independent 
of p, as expected for a particle whose mass arises 
from staggered flavor symmetry breaking effects 
of order 0(a 2 ). 



Table 2 

Results extrapolated to m q 



0. Errors are de- 



rived from 


a 250 sample 


bootstrap 




P 


6.0 


6.2 


6.4 


N/p 


1.44 ± .05 


1.37 ± .03 


1.23 ± .05 


N/U 


11.0 ± .6 


11.3 ±.4 


10.7 ± .6 


A/m 2 


.88 ± .04 


.84 ± .03 


.76 ±.07 


a _1 (m p ) 


1.91 ± .06 


2.61 ± .05 


3.39 ± .11 


a _1 (mjv) 


1.62 ±.07 


2.32 ±.06 


3.41 ± .11 




1.75 ±.06 


2.64 ± .07 


3.55 ±.16 
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